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I. Introduction 

Study of the effects of quarks and gluons confined within hadrons on 
nucleon dynamics is of fundamental importance in understanding the strong 
interaction. These effects are important for describing the short-range (SR) 
part of the A^A^ interactions. They can also be very important for describing 
the dynamics of dense nucleon matter especially at large densities. As has 
been shown in Ref. [1], the internal quark structure of baryons can have 
substantial effects on the composition and structure of neutron star matter. 
Since a nonperturbative treatment of quantum chromodynamics (QCD) is 
not possible today, the effects of the internal quark structure of hadrons on 
low-energy nucleon dynamics are investigated by using quark models (see, for 
instance, Refs.[2-9]). However, these models completely miss the effects of 
retardation in the gluon-exchange interaction, despite the fact that they can 
be important [10], and on the whole the effects of quark-gluon retardation on 
nucleon dynamics are poorly understood at present. 

The effects of quark-gluon retardation on nucleon dynamics differ 
profoundly from the well-known meson-retardation effects that are taken into 
account, for example, in the Bonn model [11]. Meson retardation effects give 
rise to nonlocality in time of the NN interaction, and hence to an energy 
dependence of the effective potentials describing these interactions, which can 
have significant effects on three- and many nucleon results [12]. Nonlocality 
in time of such interactions is an expression of a loss of probability from 
the two-nucleon subspace of the Hilbert space of hadron states. Obviously, 
quark-gluon retardation, that has to be taken into account in describing hadron 
dynamics, should also result in nonlocality in time of hadron interactions. 
However, due to confinement this must not lead to a loss of probability from 
the Hilbert space of hadron states. On the other hand, as is well known. 



of a quantum system may be nonlocal in time only in the case when the 
system is not closed, and as a result the evolution is not unitary. Similar 
problems arising due to nonlocality in time of the NN interaction has been 
discussed in [13-20]. To solve these problems quantum dynamics need to be 
extended to describe the unitary evolution of a closed system in the case 
where the interaction generating the dynamics of the system is nonlocal in 
time. For the first time, this problem has been solved in Ref.[21], where it 
has been shown that the use of the Feynman approach [22,23] to quantum 
theory in combination with the canonical approach allows such an extension 
of quantum dynamics. The generalized quantum dynamics (GQD) developed 
in this way has been shown [24,25] to open new possibilities for describing 
hadron dynamics. 

In the present paper the GQD is used for investigating quark-gluon 
retardation effects on low-energy hadron dynamics. We show that retardation 
from quark confinement results in an anomalous off-shell behavior of the 
two-nucleon amplitudes and in a lack of continuity of the evolution operator 
describing low-energy hadron dynamics that for this reason is not governed by 
the Schrodinger equation. In Sec. II we discuss some problems in describing 
hadron dynamics that arise due to quark-gluon retardation effects. The 
principal features of the GQD are reviewed in Sec. Ill, and in Sec. IV we 
show that the GQD provides an extension of Hamiltonian dynamics which 
can describe the evolution of quantum systems with confined degrees of 
freedom. The existence of such degrees of freedom gives rise to a peculiar 
dynamical situation that allows one to conclude that retardation effects of 
quark confinement on low-energy hadron can be significant. In Sec. IV this 
fact is illustrated by the example of the dynamics of nucleons with internal 
structure described by a constituent quark model. In Sec. VI we use the GQD 
to construct a model that is a generalization of the separable-potential model 
to the case when the NN interaction is nonlocal in both space and time. 
We show that despite its simplicity the model, in which retardation effects of 
quark confinement are taken into account, yields the nucleon-nucleon phase 
shifts in good agreement with experiment. Analyzing the off-shell behavior of 
the two-nucleon amplitudes in the case where the NN interaction is nonlocal 
in time, we find that quark-gluon retardation can have significant effects on 
three- and many-nucleon results. 



II. Retardation effects and low-energy hadron dynamics 



According to QCD, the physics of the strong interaction exhibits different 
relevant excitations at distinct length (or momentum) scales. At short 
space-time distances (r << 1/m) the relevant degrees of freedom are quarks 
and gluons, effectively unconfined due to asymptotic freedom. At large 
distances (r >> 1/m), on the other hand, hadronic degrees of freedom should 
be considered as relevant degrees of freedom. At the same time, due to 
confinement there are no observables, which can be associated with the quark 
and gluon degrees of freedom at least in the low-energy regime. However, in 
principle this does not mean that the quark and gluon degrees of freedom 
are not relevant in this regime. The above means only that there are 
no observables associated with quarks and gluons in the low-energy regime. 
To clarify this point, note the following. The most fundamental aspect of 
quantum theory in its present interpretation is its probabilistic character. The 
description of a physical system requires two kinds of elements, the observables 
and the states of the system. The states in turn are described by the vectors 
of a Hilbert space. The state vector determines the probabilities of finding 
some values of the observables, when a measurement is performed, and can be 
expanded in terms of the eigenvectors corresponding to the eigenvalues of a 
complete system of commuting observables. From this and the phenomenon of 
quark confinement, according to which quarks and gluons are not observable 
in free states, it follows that vectors describing low-energy states of a system 
of strongly interacting particles can be expanded in terms of eigenvectors 
associated with hadronic degrees of freedom, i.e. they belong to the Hilbert 
space of hadron states. The time evolution of the system is described by the 
evolution equation 



where \ip{t) > is a state vector of the Hilbert space of hadron states, 
and U (t, to) is the evolution operator defined on this space. Since due 
to confinement quark-gluon retardation effects must not lead to a loss of 
probability from the Hilbert space, U (t, to) must be a unitary operator 



Here we use the interaction picture. Thus due to confinement the time 
evolution of a system of strongly interacting particles at low energies can be 
described by the unitary evolution operator defined on the Hilbert space of 
hadron states. Nevertheless, the quark and gluon degrees of freedom may be 
relevant in describing low-energy hadron dynamics, since they should manifest 




(1) 



U+{t,to)U(t,to) = U(t,to)U+{t,to) = 1. 



(2) 



Let us now show that quarks and gluons confined within hadrons 
can have substantial effects on low-energy hadron dynamics. One of the 
fundamental requirement of quantum theory is that the evolution operator 
must satisfy the composition law 

U{t,t')U{t'M) = U{t,to)^ U{to,to) = 1. (3) 

In the case of an isolated system, the evolution operator in the Schrodinger 
picture Us{t2,ti) = exp{—iHQt2)U {t2,ti)exp{iHQti) depends on the difference 
(^2 — ^i) only, so that the operators V{t) = Us{t, 0) constitute a one-parameter 
group of unitary operators, with the group property 

V{h + 12) = V{ti)V{t2), V{0) = 1. (4) 

Here and below, we use the units in which c = h = 1, and Hq is the free 
Hamiltonian. If the evolution operator is assumed to be strongly continuous, 
i.e. if 

hm \\Vit2)\^>-Vih)\^>\\=0, (5) 

then from Stone's theorem it follows that this one-parameter group has a 
self- adjoint infinitesimal generator H: 

Vit) = exp{-iHt), id/dtV[t) = HV[t). 

Identifying H with the total Hamiltonian as usual, we get the time-dependent 
Schrodinger equation: i^^^^ = H\iljs{t) >, where >= V{t)\ilJs{t = 0) > 

is a state vector in the Schrodinger picture. Thus, in this case, the interaction 
being described by the interaction Hamiltonian is necessarily instantaneous. 
Obviously, the infinitesimal generator H being the operator of the total energy 
cannot contain terms depending on energy i.e. on its spectral parameter. From 
this it follows that, if the continuity condition (5) is satisfied, then the time 
evolution of a system can be unitary only in the case when the interaction 
generating the dynamics of a quantum system is instantaneous. Thus the 
evolution operator describing low energy hadron dynamics cannot be strongly 
continuous, since due to quark-gluon retardation effects hadron interactions 
should be nonlocal in time, and, on the other hand, due to confinement the 
evolution operator must be unitary. This means that quark-gluon retardation 
and confinement must result in the fact that the evolution of hadron systems 
cannot be governed by the Schrodinger equation. 

We have shown that quark-gluon retardation effects give rise to the fact 
that the evolution operator describing hadron dynamics cannot be strongly 



requirement that the evolution operator must be strongly continuous is used 
as a fundamental assumption (due to Stone's theorem it is equivalent to 
the assumption that the evolution of a quantum system is governed by the 
Schrodinger equation), this requirement is not necessary on physical grounds, 
and it is enough to require that 

< Hvit2)\^i >,^< Hviti)\^i > (6) 

for any physically realizable states > and \ip2 > [26]. Note, in this 
connection, that there are normalized vectors in the Hilbert space that 
represent the states for which the energy of a system is infinite. Such states 
cannot be considered as physically realizable, and hence the corresponding 
matrix elements of the evolution operator need not be continuous. From this 
it follows that in principle Hamiltonian dynamics can be extended to the case 
where only the physical continuity condition (6) is satisfied. As has been 
shown in Ref.[21], this really can be done by using the Feynman approach to 
quantum theory in combination with the canonical approach. As we show 
bellow, the GQD developed in this way allows one to take into account 
quark-gluon retardation in describing hadron dynamics. 

III. Generalized Quantum Dynamics 

In the GQD the following assumptions are used as basic postulates: 

(i) The physical state of a system is represented by a vector (properly 
by a ray) of a Hilbert space. 

(ii) An observable A is represented by a Hermitian hypermaximal 
operator a. The eigenvalues a^- of a give the possible values of A. An 
eigenvector > corresponding to the eigenvalue represents a state in 
which A has the value a^. If the system is in the state \ip >, the probability Pr 
of finding the value for A, when a measurement is performed, is given by 

Pr =< ^PlPvA^P >=EI <^['V> I', 

s 

where Py^ is the projection operator on the eigenmanifold Vr corresponding to 
ttr, and the sum Eg is taken over a complete orthonormal set \ipl^'> > (s=l,2,...) 
of Vr- The state of the system immediately after the observation is described 
by the vector Py^li^ > ■ 

These assumptions are the main assumptions on which quantum theory 
is founded. In the canonical formalism these postulates are used together 
with the assumption that the time evolution of a state vector is governed 



used. Instead the assumptions (i) and (ii) are used together with the following 
postulate. 

(iii) The probability of an event is the absolute square of a complex 
number called the probability amplitude. The joint probability amplitude 
of a time-ordered sequence of events is product of the separate probability 
amplitudes of each of these events. The probability amplitude of an event 
which can happen in several different ways is a sum of the probability 
amplitudes for each of these ways. 

The statements of the assumption (iii) express the well known law for 
the quantum-mechanical probabilities. Within the canonical formalism this 
law is derived as one of the consequences of the theory. However, in the 
Feynman formulation of quantum theory this law is directly derived starting 
from the analysis of the phenomenon of quantum interference, and is used 
as a basic postulate of the theory. According to the assumption (iii), the 
probability amplitude of an event which can happen in several different ways is 
a sum of contributions from each alternative way. In particular, the amplitude 

< '^2|^(^, ^o)!"^! > can be represented as a sum of contributions from all 
alternative ways of realization of the corresponding evolution process. Dividing 
these alternatives in different classes, we can then analyze such a probability 
amplitude in different ways. For example, subprocesses with definite instants 
of the beginning and end of the interaction in the system can be considered as 
such alternatives. In this way the amplitude < '02|t^(^, ^o)!^^! > can be written 
in the form [21] 

< V^2|t/(t,to)|^i >=< H^^i > + Cdt2 t'dti < V'2|5'(t2,ti)|V'i >, (7) 

Jto Jto 

where < ip2\S{t2,ti)\ipi > is the probability amplitude that if at time ti 
the system was in the state \ipi >, then the interaction in the system 
will begin at time ti and will end at time t2, and at this time the 
system will be in the state \ip2 > ■ More precisely, < ip2\S{t2,ti)\ipi > 
is the density of the probability amplitude of the above event, and 
as contributions to < 11^2^(1, to)\ipi > one has to consider amplitudes 

< HUeit2,h)\iPi >= !ll^'dt'>,!ll^'dt'^e{t'^-t[) < ^lj2\S{t'^A)\^i >, where 
e « t2 — ti and 

1, for T > 
0, for r < 0. 

The amplitude < V^2|t4(^2, ^i)|V^i > is the amplitude that if at time ti the 
system was in the state |V^i >, then the interaction in the system will begin 
in time internal ti, ti -\- e and will end in time interval t2, t2 + e, and at the 



an operator- valued generalized function of and t2 [21]. Nevertheless, it is 
convenient to call S{t2^ti) an "operator", using this word in generalized sense. 
In the case of an isolated system the operator 5(^2,^1) can be represented in 
the form [21] 

S{t2, h) = exp{iHot2)f{t2 - ti)exp{-iHoti). (8) 

To clarify the role, which the operator S{t2, h) plays in the GQD, note 
the following. The Feynman formulation is based on the assumption that the 
history of a system can be represented by some path in space-time. From the 
postulate (iii) it then follows that the probability amplitudes of any event is a 
sum of the probability amplitudes that a particle has a completely specified 
path in space-time. The contribution from a single path is postulated to 
be an exponential whose (imaginary) phase is the classical action (in units 
of h) for the path in question. In the GQD the history of a system is 
represented by the version of the time evolution of the system associated 
with completely specified instants of the beginning and end of the interaction 
in the system. Such a description of the history of a system is more 
general and requires no supplementary postulates like the above assumptions 
of the Feynman formulation. On the other hand, the probability amplitudes 
< '^2|*S'(^2, ti)\'4'i >5 in terms of which we describe quantum dynamics, are used 
in the spirit of Feynman's theory: The probability amplitude of any event is 
represented as a sum of these amplitudes. 

As has been shown in Ref.[21], for the evolution operator U{t^tQ) given 
by (7) to be unitary for any times Iq and the operator S{t2,ti) must satisfy 
the following equation: 

{t2-ti)sit2,ti) = j^yui'yhiu - h)sit2,u)sit3,ti). (9) 

Note that, since S(t2,ti) may be only an operator-valued distribution, in 
general the production S{t2,t4)S(t-i,ti) is not defined at ^4 = ^3. However, 
this does not lead to any problems because in Eq.(9) the above production 
is multiplied by the factor — t^. Eq.(9) allows one to obtain the 
operators S(t2,ti) for any ti and t2, if the operators S{t2,t[) corresponding to 
infinitesimal duration times t = t'2 — t[ of interaction are known. Thus this 
equation allows one to obtain the contributions to the evolution operator from 
the processes associated with any duration times of interaction if those from 
the processes associated with infinitesimal duration times of interactions are 
known. It is natural to assume that most of the contribution to the evolution 
operator in the limit t2 — ^ h comes from the processes associated with an 



by Hint{t2,ti), the operator S{t2,ti) can be represented in the form 

S{t2,h) = Hir^tit2:tl) + Si{t2,ti), (10) 

where Si{t2,ti) is the part of the operator S{t2,ti) which in the hmit ^2— ^^i 
gives a negligibly small contribution to the evolution operator in comparison 
with Hintih, ti)- Within the GQD the operator Hinf {t2, ti) plays the role which 
the interaction Hamiltonian plays in the ordinary formulation of quantum 
theory: It generates dynamics of a system. This operator can be regarded as a 
generalization of the interaction Hamiltonian, and it was called the generalized 
interaction operator. Obviously, the operator i^mi(^2,^i) must satisfy Eq.(9) 
in the limit t2 — ^ ti 

F(t2,ti)-0, (11) 

where 

F{t2, ti) = — (^2 — ti)Hint{t2, ti) + / / 6/^3(^4 — h)Hintit2, t^Hint{h, ti). 

If Hint{t2,ti) is specified, Eq.(9) allows one to find the operator S{t2,ti). 
Formula (7) can then be used to construct the evolution operator U {t, to) and 
accordingly the state vector 



\^{t) >= \^{to) > + fdt2 r dtiSit2,h)\^ito) > (12) 
at any time t. The corresponding evolution operator is of the form [21] 

27r J-oo [z — Ho){z — Ho) 

where z = x -\- iy, x and y are real, and y > 0, and the operator T{z) is defined 
by 

T{z)=iJ^ dTexp{izT)T{T), (14) 

with T(r) = exp{—iHot2)S{t2^ti)exp{iHoti). Thus Eq.(9) can be regarded as 
an equation of motion for states of a quantum system. It should be noted that 
in the GQD the continuity condition (6) is not used as a basic postulate, and 
plays only the role of a consistency condition of the theory: Matrix elements 
of the evolution operator (7) obtained by solving Eq.(9) must be continuous 
for any physically realizable states. It is reasonable to consider such states as 
the states for which ||ii/o|^ > 1 1 < 00. Thus the condition (6) must be satisfied 
for all vector IV' >^ ^(^o), "^(^o) being domain of Ho- 

From the mathematical point of view the requirement that Hint(t2,ti) 



means that the operator Hint{t2^ti) must have such a form that Eq.(9) has a 
unique solution having the following behavior near the point t2 = ti : 



Sit2,h)^^^ HUt2M) + o{T'), (15) 

where r = t2 — ti and the value of e depends on the form of the operator 
Hintih, ti). Since S{t2, h) and Hint{t2, ti) are only operator-valued distributions, 
the mathematical meaning of the condition (15) needs to be clarified. We will 
assume that the condition (15) means that 

< H (^[^dt2 jlytiS{t2,ti)) iV'i > ^^^^ 

< V^2| dt2 dtiHMit2. ti) j > +o(r"+2), T = t-to, (16) 

for any vectors > and \1Jj2 > of the Hilbert space. Note that the condition 
(11) has to be considered in the same sense. 

To clarify the above statement, note that Eq.(9) is equivalent to the 
following equation [21]: 

^ = T{z)G('\z)T{z), (17) 

with 



2' 

nl 



{z — Er, 

Here n stands for the entire set of discrete and continuous variables that 
characterize the system in full, and |n > are the eigenvectors of the free 
Hamiltonian: HqIti >= En\n > . Thus, instead of solving Eq.(9), one can solve 
Eq.(17) for the operator T{z). The operator T(r) and correspondingly the 
operator S{t2,ti) can then be obtained by using the Fourier transform 

i foo 

< n2|r(r)|ni >= -— / dxexp{-izT) < n2\T{z)\ni >, (18) 

27r "'-00 

where z = x + iy, x and y are real, and y > 0. At the same time, the operator 
T{z) can be directly used for constructing the evolution operator. According 
to (14) and (15), the operator T{z) has the following asymptotic behavior for 
\z\^ 00 : 

< n2\T{z)\ni > ^ < n2\B{z)\ni > +o{\z\-^), (19) 

\z\ -^00 

where 

B(z) = i I dTexp(izr)Wj(r). 



[3=1 + 6, and H-'^f{t2 - h) = exp{-iHQt2)H,,ntit2,ti)exp{iHoti) is the 
generalized interaction operator in the Schrodinger picture. From (17) and 
(19) it foUows that the operator B{z) must satisfy the following asymptotic 
condition: 

d < n2\B(z)\ni > x^f2)/ \ r^/ m n \-3\ 

dZ |2:|^C!0 

The above requirements, which the generalized interaction operator has to 
meet, mean that B{z) must be so close to the solution of equation (17) 
in the limit |z|^oo that this differential equation has a unique solution 
having the asymptotic behavior (19). The operator B[z) represents the 
contribution which H^^]{t) gives to the operator T{z), and was called the 
effective interaction operator. It should be also noted, that the operator T{z) 
satisfies the equation 

rpf \ \ ( ^^ T{z2)\n><n\T{zi) 

T{zi) - T{z2) = {Z2 - zi) E — — — — -, (20) 

n [Z2 — rjn)\Zi — tjn) 

provided that S{t2,ti) satisfies Eq.(9). 

As has been shown in Ref.[21], the dynamics governed by Eq.(9) 
is equivalent to Hamiltonian dynamics in the case where the generalized 
interaction operator is of the form 

Hint{t2M) = -2iS{t2-h)Hi{ti), (21) 

Hiiti) being the interaction Hamiltonian in the interaction picture. In this 
case the state vector |'0(t) > given by (12) satisfies the Schrodinger equation 

The delta function 5{t) in (21) emphasizes that in this case the fundamental 
interaction is instantaneous. Thus the Schrodinger equation results from the 
generalized equation of motion (9) in the case where the interaction generating 
the dynamics of a quantum system is instantaneous. At the same time, Eq.(9) 
permits the generalization to the case where the operator Hjnt{t2-,'ti) has no 
such a singularity as the delta function at the point t2 = ti [21]. In this case 
the fundamental interaction generating the dynamics of a quantum system is 
nonlocal in time: The evolution operator is defined by Hint{t2^ ti) as a function 
of the time duration r = t2 — ti of the interaction. In a more general case, the 
generalized interaction operator has the following form [25]: 



where the first term on the right-hand side of (22) describes the instantaneous 
part of the interaction generating the dynamics of a quantum system, while 
the term Hnon represents its nonlocal-in-time part. In order that the dynamics 
governed by Eq. (9) be different from the dynamics governed by the Schrodinger 
equation with the interaction Hamiltonian Hj{t)^ the nonlocal operator Hnon 
must satisfy the condition 

{t2 — ti)Hnon{t2iti) ^ I dt4 dt^[t4^ — t^)Hnon{t2it4)Hnon{t?,itl) i ^2" ^^l- 

Jtl Jtl 

(23) 

In fact, in the local case, where Hint{t2,ti) is of the form (21), the first term in 
the operator F{t2,ti) is zero, and from the condition (11) it follows that the 
second term in this operator must tend to zero in the limit t2 ti so rapidly 
as it is required for Eq.(9) to have a unique solution. In the nonlocal case, 
both terms in the operator F{t2-,ti) are not zero, and, as we show in Sec. IV, 
they do not decrease when t2 ti so rapidly as it is needed. Only their 
sum can vanish fast enough for the condition (11) to be satisfied. From this 
it follows that the nonlocal part of the generalized interaction operator must 
satisfy the condition (23). 

The principal feature of the GQD is that it provides an extension of 
Hamiltonian dynamics which can describe the unitary evolution of a quantum 
system with confined degrees of freedom. In the next section we will show that 
the existence of such degrees of freedom, that can manifest themselves only 
in some retardation of the interaction in the system, gives rise to a peculiar 
dynamical situation that allows one to conclude that retardation effects of 
quark confinement on low-energy hadron dynamics can be very significant. 

IV. Effects of confined degrees of freedom on the dynamics 
of a quantum system 

Let us consider the evolution problem for two nonrelativistic particles in the 
c.m.s. We denote the relative momentum by p and the reduced mass by ji. 
Assume that the generalized interaction operator in the Schrodinger picture 
h\^^[t) has the form 

< P2|//S(r)|pi >= ^\v2Mvi)f{r), (24) 

where /(r) is some function of r, and the form factor (/?(p) has the following 
asymptotic behavior for |p| — ^ oo : 



Let, for example, f{p) be of the form 

^(P) = T^ + ^(P), (26) 

IpI 

and in the hmit |p|^oo the function g{p) satisfies the estimate g{p) = 
o(|p| '^), where 5 > a, (5 > |. In this case, the problem can be easily 
solved by using Eq.(15). Representing < p2|r(z)|pi > in the form 
< p2|T(z)|pi >= ip(jp2)(p* ipi)t{z) , from (15), we get the equation 

with the asymptotic condition 

t(z),,^ h{z) + o{\z\-% (28) 



-OO 



J 2 

fi{z) = i Jq° drexp{izT)f{r), and Ek = The solution of Eq.(27) with the 
initial condition t{a) = Qa, where a G (— oo,0), is 

tiz) = <,„(! + (.- a),. / A (^_gf_^^) ) ' . (29) 

In the case a > |, the function t{z) tends to a constant as z— ^ — oo: 

t{z) X. (30) 



-oo 



Thus in this case the function fi{z) must tend to A as z^—oo. From this it 
follows that the only possible form of the function /(r) is 

/(r) = -2iA(5(r)+/'(r), 

where the function /'(r) has no such a singularity at the point r = as the 
delta function. In this case the generalized interaction operator h\^^[t) has 
the form (21) and hence the dynamics generated by this operator is equivalent 
to the dynamics governed by the Schrodinger equation with the separable 
potential 

< P2|i^j|Pl >= A(^*(p2)(^(pi). (31) 

Solving Eq.(27) with the boundary condition (30), we easily get the well-known 
expression for the T matrix in the separable-potential model 



< P2|r(z)|pi >= A(p*(p2)^(pi) fl - A / d^A^M^) 



(32) 



Ordinary quantum mechanics does not permit the extension of the above 
model to the case a < ^. Indeed, in the case of such a large-momentum 
behavior of the form factors v?(p), the use of the interaction Hamiltonian given 
by (31) leads to the ultraviolet divergences, i.e. the integral in (32) is not 
convergent. We will now show that our formalism allows one to extend this 
model to the case < a < |. Let us determine the class of the functions fi{z) 
and correspondingly the value of (3 for which Eq.(27) has a unique solution 
having the asymptotic behavior (28). In the case a < ^, the function t{z) 
given by (29) has the following behavior for —oo : 



tiz) ^ -^_J,{-zr-2 + 62(-^)'""' + o(|z|2"-i), (33) 



where bi = -|cos(a7r)7r ^|ci| ^(2//)" 2 and 62 = ^il^l^ " - bl{M{a) + g^^) 
with 

12 Icip 



l^(k)|^- 



M(a)=/^J^A 



a- Ek 

The parameter 61 does not depend on Qa- This means that all solutions 
of Eq.(27) have the same leading term in (33), and only the second term 
distinguishes the different solutions of this equation. Thus, in order to obtain 
a unique solution of Eq.(27), we must specify the first two terms in the 
asymptotic behavior of t{z) for z — ^ —00. From this it follows that the functions 
fi{z) must be of the form 

hiz) = 6i(-^)"-^ + 62(-^)'"-\ (34) 

and (3 = 2a — 1. Correspondingly the functions /(r) must be of the form 

/(r) = air-"-^ + a2r-2^ (35) 

with ai = —ibir^^{l — 2a)exp[i{—^-^ + j)7T], and 02 = —b2r^^{l — 2a)exp{—ia7T), 
where T'(z) is the gamma-function. This means that in the case o. < ^ the 
generalized interaction operator must be of the form 

< p2|i?S(r)|pi >= /(r)(^*(p2)(/p(pi) = 

= ai^*(p2)^*(pi)r— ^ + a2(^*(p2)^(pi)r-2", (36) 
and, as it follows from (29) and (33), for the T matrix we have 

< P2|T(^)|pi >= 7V(z)^*(p2)^(pi), (37) 

with 



where Qa = h\ ^6i|a|5~" + a2r(l — 2a)exp{ia'K) — h\M{a)^ . It can be easily 
checked that N{z) can be represented in the following form 

N{z) = ^ -. 

-b2 + hi- z)2-^ + M{z)bl 

By using (13) and (37), we can construct the evolution operator 
< P2|t/(i,io)|Pi >=< P2IP1 > -7:;- / dx 

[z — rjp^){z - tjp^) 

where z = x + iy, and y > 0. The evolution operator U{t^to) defined by (38) 
is a unitary operator satisfying the composition law (3), provided that the 
parameter 62 is real. 

We have stated the correspondence between the form of the generalized 
interaction operator and the large-momentum behavior of the form factor 
i^ip). In the case a > ^, the operator H^^^(t) would necessarily have the 
form (21). In this case the fundamental interaction is instantaneous. In 
the case < a < ^ (the restriction q; > is necessary for the integral in 

(29) to be convergent), the only possible form of i7f^|(r) is (36), and hence 
the interaction generating the dynamics of the system is nonlocal in time. 
Thus the interaction generating the dynamics can be nonlocal in time only 
if the form factors have the "bad" large-momentum behavior that within 
Hamiltonian dynamics gives rise to the ultraviolet divergences. 

Let us now show that in the case ot < \ the dynamics generated by 
the generalized interaction operator (36) is not equivalent to Hamiltonian 
dynamics. The evolution operator U (t, ^o) given by (38) satisfies the 
composition law (3) for any ^q, t and t' ^ provided the operator S{t2,ti) satisfies 
Eq.(9) for any ti and ^2- In the case of model under study the operator 
S{t2,ti) is of the form 

Sit2,ti) = D{t2,h)F{t2,h), (39) 

where F{r) is a function of r, and D{t2,ti) = exp{iHQt2)\(p >< (p\exp{—iHQti) 
is the operator-valued distribution such that 

< P2|^(^2,ii)|pi >= exp{iEp^t2) < P2|^ >< V'Ipi > exp{-iEp^ti) , (40) 

with < (p\p >= (fijp). By using (39) and (40), for the function -F(r) we get the 
following equation: 
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X < (p\k >< k\(p > F{t2 - U)F{U - h). (41) 

Let us examine this equation in the Umit ^2 — ^ • For this we have to change 
the variables: Ti = 6i/u, k = z/q,y, i = 1,2,3,4. With such a change of the 
variables, Eq.(41) can be rewritten in the form 

{e2-ei)Fie2/iy^-ei/iy') = f^'^ de, fUe^{e^-e^) j d'kexp[-iE,^ie,-es)]x 

X < ^|q, X q,|^ > F{02/i^^ - e^/v'^)F{e^/v^ - Oi/v^), (42) 

3 _ 

Iqi, > being the basis vectors such that < k|qi, >= v~^5{k — v^p)- The 
completeness condition for these basis states reads 

j d^qy\i\y >< = 1. (43) 

Thus |k > and |qj, > are eigenstates of the operator Hq corresponding to 
distinct momentum scales. Since < (p\c{v >= ci|z/q,y|~" + o(z/~2~") when 
ly^oo, by letting t2^ti, from (42) we get 

{02 - 0i)F{02/iy' - Oi/u^) = r de^ d0sFi02/iy' - O^/u^) x 

X F{e,/u^ - ei/u^)I{0^ - ^3, z/), z/ ^ oo, (44) 

with 

1(0, u) = \ci\^e I (iVa:p(-i^,.^)|q.|"'" + o(z/-i+"). 

This expression emphasizes that the relevant momentum scale of intermediate 
states in Eq. (42) tends to infinity as ^2 — ^ • In order to demonstrate that 
Eq.(44) plays the key role in the case a < ^, note that from this equation it 
follows that in the limit t^O the function F[t) behaves like aiT~^~^ with 

ai = {02 - ^i)^-" de, /^'' de,i{e, - ^3, 00)) ' = 

= ^7r->ir2(2//)"-ir-i(l - 2a) cos{a7r) exp[i{-^ + hn], (45) 

i.e. this equation determines the main term of the generalized interaction 
operator (see Eq.(36)). We see from (45) that only the intermediate states 
with infinite momentum contribute to the parameter ai. This manifests itself 
in the fact that this parameter depends only on the leader term in the form 
factor (25). Note in this connection that for any vector \ip > of the Hilbert 
space H represented in the form 



with V^(k) =< k.\ilj >, one can construct another vector 

>= J ^(k/z/)z/-i|k > d^k, (47) 

that represents the same physical state, if we scale k — ^ z/k, i.e. 

IV^^ >= J V^(qi.)|qi. > (i^g^., (48) 

where = k/z/. Varying the parameter z/, we get a set of vectors having the 
same norm 

> W = {J rMi^Md'qS'- (49) 

Each of the vectors li/j^y > belong to the Hilbert space Ti even when the 
parameter u is letting to infinity. To clarify this point, let us consider the 
state IV^ >= / ?/^feg(k)|k > d^k, where ^fco(k) is zero everywhere outside the 
subset A(A;o) (£^fco ^ < E^q + cEk^, e <^ 1) of the spectrum of Hq. Such 
a state is eigenstate of the projection operator PA(fco) subset A{ko). 

Correspondingly the state 

IV^^ >= J tlJq,iq^)\q^ > d^qy (50) 

is eigenstate of the projection operator on the subset A (go) {Eq^ < E^ < 
Eq^ + eEq^^, Eq^ = Ek^u~'^). The projection operators are defined for any 
subsets even when their location tend to the infinite part of the spectrum. 
From this it follows that even when v tends to infinity eigenvectors of the 
projection operator -PA(go) belong to the Hilbert space H. However, for 
describing such states we have to make a change the scale by letting u to 
infinity. For any u the vectors [if^y > can be represented as vectors of the 
Hilbert space L?{My) of square integrable functions '^(Qz^)? where denotes 
the momentum-space corresponding to the scale u. Obviously, for any finite u 
the spaces L?{Mi) coincides each with other, since '0(Qi^) = V'X^) — '0(k/^)- 
This, of course, is not true for the space L'^{Moo) of square integrable functions 
'^{q.in), with being the momentum in the case when the scale tends to 
infinity. Since for any |k| < cxd 

JWL < ^^^^ I ^(^'^) < ^l^'^ > ^^^'^ = 

= JWL j ^(^'^) < > i^~~'d^qu = \mi^^ik/iy)p--^ = 0, (51) 

each of the vectors \ipin > of L^(M(x,) is orthogonal to all vectors of L^(Mi). 
Thus the Hilbert space describing the states of the system under study is 



where Tip is the space that can be reahzed as the space L^(Mi), and TYm is 
the space that can be reahzed as the space L2(Moc). The manifold of the 
physically realizable states, being the domain of iifg, is dense in Tip. On the 
other hand, the vectors of Tiin represent the states with infinite energy, and 
are not physically realizable. 

Let us now show that the matrix element < Ojl^/^^^ >, where 

R{t^ to) is defined by U (t, ^o) = 1 + iR{t^ to), are not continuous at t = 0. From 
(7) and (8) it follows that the operator R{t, to) can be represented in the form 



R{t,to) = —i f dt2 [ dtiexp{iHot2)T(t2 — ti)exp{—iHoti). (53) 

By using (53), we can write 

< V'?) \R{t, 0) I V'?) >= -i dt2 J^' dti I d^k' j d^kexp{iEk't2) x 

X exp{-iEkti) < ^l^^\k' >< k'\f{t2 - ti)\k >< klV'W > . (54) 
Taking into account (25), (40) and (47), and letting u to infinity, we get 

< 4l^\R{t,0)\4n >= JWL < ^l^^m,0)\4'^ >= -^«iki|%lim /J''c?^2X 
X J^'dOi J d^ql J d^q^exp{iEq,^02)exp{-iEqJi)ilj;{c^^)ilji{q^)x 

xz/ >q^| "|z/qj " I 1 

where 6i = tiu'^, ipi{ci,^) =< cii/\ipl >, i = 1,2. From this it follows that in the 
limit z^— ^oo the matrix elements < iljl^'>\R{t,0)\ipl^'> > are scale invariant, and 
we have 

< 4n\Rit,0)\4i^ >= -iai\ci\^ de2 fjdOi j d'q'J d^qexp[iE,,e2)x 

X exi?(-zE,^i)V^2*(qOV^i(q)|qnqr"(^2 - ^i)"^"". (55) 

Thus the matrix elements < ilji'^\R{t, ^)\i^^ > are independent of and do not 
tend to zero as t — ^ 0. This means that the amplitudes < '0-^''|i?(t, 0)|?/^£'' > 
are not continuous at t = since, as it follows from the definition, R(0, 0) = 0. 
From this it follows that the evolution operator (38) is not weakly continuous. 
On the other hand one can show that for the physically realizable states 
Itjji > and |'02 > the matrix elements < ?/^2|f^(^, O)]?/'! > tend to < '02|V^i > 
as t — ^ and hence the condition (6) is not violated. Nevertheless, the 
evolution operator V(t) = Us(t,0) is not continuous and hence the group of 



that in this case the time evolution of a state vector is not governed by the 
Schrodinger equation. The cause of this discontinuity is quite obvious. In 
fact, from the point of view of the states with infinite energy any time interval 
6t is infinite and hence the matrix element < 0)|?/^£^ > must be 

independent of 6t^ i.e. must be constant. 

The fact that in principle for describing the evolution of a quantum 
system during infinitesimal time intervals one has to take into account the 
intermediate states with infinite energy is a consequence of the principle of 
uncertainty. In the case of Hamiltonian dynamics, nevertheless, such states 
give no contributions because the matrix elements < p2|t/(t2, ^i)|pi > vanish 
sufficiently fast when momenta tend to infinity. In this case the space Tim can 
be ignored, and the dynamics can be restricted to Tip. However, in general 
case the matrix elements < p2\U{t2^ti)\pi > may have such a large-momentum 
behavior that one cannot ignore the space Tim in describing the dynamics of a 
quantum system. To better illustrate this point, let us consider the operator 
Vp{t) = PV{t)P, P being the projection operator on the subspace Tip. Assume 
that the group of the operators Vp{t) has a self-adjoint infinitesimal generator 
A. Then for |?/; >e V{A) we have 

^«f^-„-»^l<^>. (56) 

From this it follows that A = Hq + Ai, with Ai = — lim {j:exp{—iHQt)R{t^ 0)) . 
By using (40) and (53), for >G T>[Hq) and |k| < oo, we get 



<k|Ai|V^ >=ilini(^-^ dt2 dti j (fkiexp[iEk{t2 - t)\x 

xexp{-iEj,,h)if\V.)ip{\^^) < k|f (^2 - ti)|ki > V'(ki)) = 

= iai|ci|2C^*(k)lim Q/Jdi2 /J'^iii(i2-^i)-^"") =0, (57) 

where C = J (i^/ci(/?(ki)'?/^(ki). Hence A = Hq, i.e. the infinitesimal generator of 
the group of the operators Vp{t) is equal to the free Hamiltonian. This means 
that the dynamics of the system cannot be restricted to the space Hp. 

It is important to now whether it is possible to find the system at time 
t in a state G Tim, if initially at time to the state of the system was 
physically realizable. The amplitude of this possibility is given by the matrix 
element < ?/^i„|[/(t, 0)|V^i > with >£ Tip. For this matrix element we can 
write 



= aicl lim^ Jo^ J ^^^^ / d^k'^xp{^Eq^92) x 

X {02 - ^i)-^-"|q,|->*(q,)V'i(k>-i = 0. (58) 

Thus, if at time the state of the system was physically realizable, the 
probability to finding the system at time t in any state \ip >£ Tim is equal 
to zero, despite the fact that one cannot ignore the subspace Tiin- This 
means that the states |V' are not observable. This is not at variance 

with the fact that one have to take into account the space Tiin in describing 
the time evolution of the system. Indeed, in the limit all matrix 

elements < V^2|^(^, 0)|V^i > tend to zero, provided that \ipi >G Tip. However, 
for |'02 >^ Tip the matrix elements < '^2\Rit:^)\i>i > vanishes faster than 
< ilJin\R{t,0)\ilJi > when t— ^0. This results in the fact that, as we have seen, 
the intermediate states belonging to Tiin are responsible for validity of the 
composition law (3) for infinitesimal time intervals. Thus these states play the 
key role in the time evolution of the system in the infinitesimal neighborhood 
of the point t = 0. 

To clarify this point, note the following. The concept of nonlocal- in-time 
potentials were first introduced within the optical-potential model. The optical 
potentials are introduced in the case when only state vectors belonging some 
subspace of the Hilbert space A are included explicitly in the description 
of the time evolution of a quantum system. Such potentials which globally 
accounts for the coupling between the subspace A and its complementary part 
B are nonlocal in time, and depend on the history of the system. The nonlocal 
form of the optical potentials is an expression of the loss of probability from 
the subspace A, i.e. of the fact that the evolution operator defined on A is 
not unitary. The nonlocality in time of effective interaction operator in the 
case, when the dynamics of a quantum system is restricted to some subspace 
of the Hilbert space, is a consequence of the coupling between this subspace 
and its complementary part. This clarifies the fact that in our model the 
interaction operator becomes nonlocal in time in the case a < ^ when the 
ultraviolet behavior of the form factors is "bad", and one cannot ignore the 
subspace Tiin of states with infinite energy. In this case we also have to 
deal with two subspaces: Tip and its complementary part Tiin- However, in 
our model the coupling between these subspaces does not lead to the loss 
of probability from the subspace Tip : The evolution operator PU{t,to)P, P 
being the orthogonal projection on Tip, is unitary. The remarkable feature 
of the GQD is that it allows one to restrict the description of the unitary 



account the couphng between this space and Tim- In fact, as we have seen, 
the space Tim manifests itself only in the behavior of the evolution operator 
U{t,to) in the infinitesimal neighborhood of the point t = ^o- On the other 
hand, the physical meaning of the generalized interaction operator Hint{t2,ti) 
is that it determines the evolution of the system in this neighborhood. In 
order that Hint{t2, ti) contain all needed dynamical information, in the limit 
h^h the operator to) = It^dh It^ ^tiHintih^ti) must be close enough 

to the true evolution operator U{t^tQ). Thus the solution of the evolution 
problem is divided in two stages: The first step is to determine the operator 
Hint{t2-,ti) describing the time evolution in the infinitesimal neighborhood of 
the point t = to, where the states belonging to Tim play an important role, and 
the next step is to describe the dynamics restricted to Tip starting with the 
above generalized interaction operator that takes into account the coupling 
between the spaces Tip and Tiin- 

To illustrate this point, let us come back to our model. From (29) it 
follows that the solution of Eq.(27) can be represented in the form 

t{z) = ^hm^ga{l + gaJ{z,a)y^ , (59) 

where 

In the case a < |, J{z,a) tends to infinity like |a|2~" when a— ^— oo. From 
this it follows that, for the solution t{z) to be nonzero, Qa must decrease like 
|a|"~2 as a— ^ — oo. Thus the leader term of Qa in the limit a— ^ — oo must be 
of the form 6i|a|~2+"j where the parameter 61 is completely determined by 
the above requirement of nontriviality of the solution. Indeed, by using the 
identity J{zi, Z2) = J{zi, z) + J{z, Z2), we can rewrite (59) in the form 

t{z) = Jjm^Qa (1 + 9aJ{0, a) + gaJ{z, 0))"' . (60) 

From this formula we see that QaJi^^, o) must tend to one as —00, since in 
this limit Qa tends to zero. Hence we have 

61 = - lim ( |ao|^+V^- /^'^.TT^^^^r' = 

= - (|ao|^+1ci|2/^iVlV^(qzn)|')"' , (61) 

1 —- 

where a = aoz/^, ao G (-00, 0), V'(qm) = hm\~''iEq.^ - ao)~2EgJ. Thus the 



LP'{Moo) being the realization of the space Tim- From (60) it also follows that 
the next order term in the asymptotic behavior of Qa must be of the form 
62|a|~"^'^^", where the freedom of choice of the parameter 62 is bounded only by 
the requirement of unitarity of the resulting evolution operator. Thus, in the 
limit —00 the function t{a) = ga behaves like 6i|a|^2+" + 62|ci|~^^^"5 and 
hence the generalized interaction operator must be of the form (36), and the 
corresponding solution of Eq.(17) is as follows: 

< P2|T(^)|pi >= -v'(p2)f(Pi) {b2bf + 1 A^^j^g^) ' . (62) 

It is easy to verify that (62) is equivalent to the above obtained expression 
(32). Here we are dealing only with the functions belonging to the space 
L^(Mi). Thus after specifying the form of the generalized interaction operator 
the description of the dynamics can be restricted to the space of the physically 
realizable states. Note that in the case ct > ^ such problems do not arise. 
Indeed, in this case Qa in (29) must tend to A being an arbitrary real constant. 
Thus, in this case we directly get the expression (32), which is written in 
terms of functions belonging to the space L^(Mi). 

By the example of the exactly solvable model we have shown that 
nonlocality in time of the interaction generating the unitary dynamics of a 
quantum system should be associated with the existence of the space Hm of 
unobservable states that has to be taken into account in describing the time 
evolution of the system. In our model the space TY^ may be considered as 
the space of two-nucleon states. The quark structure of the nucleons being 
in such states does not directly manifest itself since, if the system is in the 
state \ip >G Hp, the probability of finding free quarks, when a measurement is 
performed, is equal to zero. In this case, the space Hin may be considered as 
the space of states in which the quark structure of nucleons manifests itself. 
The spaces Hp and Hin have the same structure. In the general case the 
structure of these spaces may be different. For example, in the case of hadron 
dynamics the space Hin of states corresponding to infinite momentum scale, 
i.e. to infinitesimal length scale, may describe any states of quarks and gluons, 
while due to confinement the space Hp describes only observable hadron states 
(color singlets). 

Note in this connection that, as we have stated, there are two different 
relevant scales in the physics of the strong interaction, and their separation is 
so significant that the length scale of processes in which quarks and gluons 
are relevant degrees of freedom can be regarded as infinitesimal compared to 
the relevant length scale of the low-energy hadron physics. Indeed, these 



intervals. From this in turn it follows that the relevant momentum scale of 
states in which quarks and gluons can manifest themselves must be infinitely 
large. Nevertheless these degrees of freedom can have significant effects on 
hadron dynamics, and for its describing one has to take into account the 
coupling between the space of hadron states and the space describing states 
of quark and gluons, and, as usual, this coupling should result in nonlocality 
in time of hadron interactions. However, in contrast with ordinary open 
quantum systems, due to confinement this nonlocality must not lead to a 
loss of probability from the hadron system: The probability of finding the 
system in a state containing quarks and gluons, when a measurement is 
performed, is equal to zero. Nevertheless one cannot ignore the subspace 
describing such states. As we have seen, it is just the case when one has to 
go beyond Hamiltonian dynamics. The remarkable feature of the GQD is that 
it provides the extension of quantum dynamics to this case. This opens new 
possibilities to find the link between the quark-gluon dynamics and low-energy 
hadron dynamics. We can try to construct the operator of hadron interactions 
by using the quark-gluon dynamics, and then use it for describing hadron 
dynamics. As we will show in the next section, for constructing this operator 
one can also use some quark models. 

The essential lesson we have learned from the previous analysis is that 
the existence of the quark and gluon degrees of freedom confined within 
hadrons gives rise to the peculiar dynamical situation: In the large-momentum 
limit the matrix elements of the evolution operator do not decrease so 
rapidly as it is required by ordinary quantum mechanics, and this results 
in the above-mentioned lack of continuity of the evolution operator, and 
correspondingly in the fact that hadron dynamics is not governed by the 
Schrodinger equation. More precisely, the above are retardation effects from 
quark confinement on hadron dynamics, since the existence of some external 
degrees of freedom in itself cannot have such effects on the dynamics of 
a quantum system. As we have seen, the retardation effects from quark 
confinement can be very essential. For example, as we will show in Sec. VI, 
these effects give rise to an anomalous off-shell behavior of the two-nucleon 
amplitudes. 

V. Generalized Interaction Operator and Quark Models 

As we have shown, quark-gluon retardation must result in nonlocality in time 
of hadron interactions. From this it follows that low-energy hadron dynamics 



generalized interaction operator describing the history of the hadron system 
in an infinitesimal time interval. Here we mean the dynamical system whose 
states correspond to a complete set of hadronic observables. In the case 
of low-energy hadron dynamics, we can consider as infinitesimal such time 
intervals, for which only the processes being the fundamental at the hadronic 
level, for example, the processes being described by the baryon-baryon-meson 
and four-baryon vertices, are relevant. The generalized interaction operator 
Hint{t2-, ti) represents the history of the hadron system in infinitesimal time 
intervals, during which the quark and gluon degrees of freedom can manifest 
themselves, and in principle it should be extracted from QCD. However, 
a nonperturbative treatment of QCD is not possible until now, and for 
constructing the generalized interaction operator we have to restrict ourselves 
to using some quark models. 

Let us now show that the dynamics of nucleons with the internal 
structure described by a constituent quark model is governed by the equation 
of motion (9) with a nonlocal-in-time operator Hint{t2-,ti)- In constituent 
quark models nonstrange and strange baryons are considered to be clusters 
made of three valence quarks. The Hamiltonian used in these models has the 
following form: 

2 

H = Y.^- Tcm. + V^onir,,) + V,,(r,,), (63) 

where r^j represents the interquark distance, Tcm. is the kinetic energy of the 
center of mass motion, Vqqirij) is the one-gluon exchange (OGE) potential, 
and Vcon(jij) is the confinement potential 

Here ac is the confinement strength, and are the SU(3) color matrices. 
In the models baryon states are represented by the wave functions of a 
three-quark oscillator with the Hamiltonian (63). As is well known, such 
models provide a very satisfactory description of the baryon spectra. In order 
to take into account the long- and intermediate- range part of the interaction, 
in some models in addition to the OGE interaction, the quarks belonging 
to different three-quark clusters interact via scalar and pseudoscalar meson 
exchange. 

For solving the two-nucleon problem within a constituent quark model, 
one has to consider the dynamics of the six-quark system. Due to confinement 
potential the six-quark system can be only in the two-baryon states lip; 2B > 
or in the baglike states \ib: 6q >. However, for large time intervals the system 



Tc, during which the system can be in baghke states, is very small and depends 
on the confinement length Ac- For the time intervals large compared with Tc, 
the matrix elements < i{j2'i6q\U{t^to)\ilJi > describing the probability to find 
the system in the state |?/^2; 6g > at time t are neglectable, and we can write 

< V'2;6g|t/(t,to)|V'i >=0. (64) 

From this it follows that for large time intervals the operator t/2s(^,^o) 
obtained by projecting t/2B(^,^o) = P2BU{t^to)P2B^ where P2B is the projection 
operator on the subspace 7^2^, can be considered to be unitary. 

Let us now show that the representation (7) is valid for this operator. By 
using the composition law (3) for the evolution operator [/(t,to), we can write 

N 

U{t,to)= I[U{t„t,_i), (65) 
i=i 

where tj = to + j^, ^ = {t — to)/N, = t. Within Hamiltonian dynamics the 
evolution operator U{tj+i,tj) has the form 

U{tj, tj-i) = exp{iHotj)exp[—iH{tj — tj-i)]exp{—iHotj-i). (66) 

In the limit ^ we can write 

U{tj,tj-i) ^ exp{iHotj+i){l — i^H)exp{—iHotj-i). (67) 

By using (65), the evolution operator can be represented in the form 

N 

Uit,tQ) = lim n exp{iHotj){l - i^H)exp{-iHQtj-i). (68) 
It is easy to verify that (68) can be rewritten in the form 

U{t, to) = 1 - Jim E i^^HIitl) + E eHi{ti)U{ti, tk)Hi{tk)\ . (69) 

^^^1=1 V k=l J 

Taking the limit ^^0, from (69), we get 
U{t, to) = 1 - r dt2 t dti {2i5{t2 - h)Hi{ti) + Hi{t2)U{t2, h)Hi{t{)) . (70) 

Thus, within Hamiltonian dynamics, the evolution operator can be represented 
in the form (7) with the following operator S{t2,ti): 

Sit2,h) = -2iS{t2 - h)Hi{ti) - Hi{t2)U{t2M)Hi{ti). (71) 

Correspondingly, for the operator ?72s(t,to) defined on the two-baryon space 
7^2S we can write 



+ / dt2 rdti<ij2,/3\S2Bit2,ti)\i;i,a>, (72) 

where 

< i^2, P\S2Bit2,ti)\i/ji,a > = < ^2, (^\P2BS{t2M)P2B\ll^U(^ > = 

= -< ^2.(^\P2Bexp{iHft2){2i5{t2 - 

+ HfUsit2, ti)Hi''^}expi-iHi%)P2B\^u a > . (73) 

Here we use the interaction picture in which the interaction in the baryon 
clusters is included into the "free" Hamiltonian Hq^'' of the corresponding 
two-baryon channel H = hI^^ + and \'ip^ot > is the vector belonging to 

the channel subspaces Tia- 

We have shown that the operator t/2B(^,^o) can be represented in the 
form (7), where the operator S{t2,ti) is given by (73). As has been shown in 
Ref. [21], for the operator t/2s(^,^o) having the form (7) to be unitary for any 
times t and to? the operator 52^(^2, ^i) must satisfy the equation 

(^2 - ti)S2B{t2, h) = dU t dh{U - h)S2B{t2, U)S2B{h, h) . (74) 

However, the operator U2B{t,to) can be considered as a unitary operator only 
for time intervals t — to large compared with Tc- Correspondingly, the operator 
S2B{t2,ti) satisfies the equation (74) only for time intervals t2 — h » Tc- 
On the other hand, due to confinement time intervals, which are relevant for 
the description of low-energy nucleon dynamics, are large compared with Tc 
(here we are dealing with two separated length scales). From this it follows 
that the low-energy dynamics of the two-nucleon system predetermined by a 
constituent quark model is governed by Eq.(9). 

Up to now we restricted ourselves to the consideration of the two-nucleon 
dynamics. On the other hand, there are time intervals large, compared with Tc, 
that are much smaller than the time intervals relevant for the many-nucleon 
interactions. Hence, for such "infinitesimal" intervals of time the operator 
SN{t2,ti), which is the two-nucleon operator on the many-nucleon Hilbert 
space, satisfies Eq.(9). For this reason, the operator SN{t2,ti) related to 
the operator S2B{t2,ti) in the ordinary way can be used as the generalized 
interaction operator generating the dynamics of the many-nucleon system 

Hint{t2M) = SN{t2M)- (75) 

From (73) it follows that this operator can be represented in the form 
(22), where the long- and intermediate-range parts of HAt) present the 



in the six-quark Hamiltonian, and its short-range part is a potential arising 
due to the one-gluon-exchange potential in this Hamiltonian. The nonlocal 
term Hnonih^ti) describes the contributions from processes in which during 
some time intervals retardation from quark confinement can take place. 
The contribution from such processes cannot be neglected, and the effective 
NN interaction extracted from a constituent quark model must contain the 
nonlocal-in-time term leading to the above-mentioned dynamical situation. 
This is an expression of the fact that there are some effects, related to 
the short-range NN system, that cannot be described in terms of the A^A^ 
potentials (see, for example, [8]). 



VI.Nonlocality in time of the NN interaction and an 
anomalous off-shell behavior of two-nucleon amplitudes. 

The results given in Sec.IV can be used for constructing models which allows 
one to take into account the quark-gluon retardation effects in describing the 
NN interaction. As the first step in this direction, we can generalize the 
ordinary separable-potential model that is widely used in nuclear physics. 
At the present time, for two-nucleon separable potentials it is usually used 
the Yamaguchi and Tabakin form factors [27,28] which in spin-triplet and 
spin-singlet channels are of the form 

grip) = 2 ! (''6) 

(p2 + 77)(p2 + ^2)fc' 

where grip) and grip) are respectively Yamaguchi and Tabakin form factors, 
and p = |p|. Such form factors with the parameters given in Ref. [29] have 
been recently used by Rupp and Tjon [30] and Adnikari and Tomio [31]. In 
Refs.[31], for example, the Tabakin form factor has been used in one of the 
nucleon-nucleon spin channels and Yamaguchi in the other. 

As it follows from the analysis of Sec.IV, the separable potentials 
cannot take into account the quark-gluon retardation in describing the NN 
interaction, since in this case the NN interaction is instantaneous. For the 
interaction to be nonlocal in time, in the large-momentum limit the form 
factor (f{p) must behave like (25) with < a < ^, i.e. must have the 
asymptotic behavior that in the case of the separable potentials leads to the 
ultraviolet divergences. In this case the generalized interaction operator must 



the coupling between the space Tip of the observable states of the system 
under study and the space Ti^ of unobservable states manifests itself only in 
the asymptotic behavior of the form factor v^(p) in the limit |p|^cxD, and 
in the behavior of the H^^^(t) in the infinitesimal neighborhood of the point 
r = 0. Note that in the case of the NN interaction the space Hp should 
be considered as the space Hn of the two-nucleon states, and Tim as the 
space Tiq of states, in which the quark structure of nucleons manifests itself. 
Since the parameter ai in (36) is completely determined by the parameters 
a and ci characterizing the asymptotic behavior of the form factor (f{p), we 
have only three free parameters a, ci, and 02 responsible for the coupling 
between the spaces T^tv and Tiq. However, this true only in the case of the 
rank-one separable approximation. In this case these three parameters can be 
determined, for example, by fitting the NN phase shifts. In a more general 
case, when we cannot restrict ourselves to this approximation, we have much 
more parameters responsible for the coupling between the spaces Hn and Tiq, 
and these parameters should be extracted from quark-gluon dynamics. Let us 
consider the rank-one separable model of the A^A^ interaction. Let the form 
factor <^(p) in (24) be of the form 

^(p) = X(P) + C2py(p), (78) 

with 

where d, ci and C2 are some constants. Since a < |, the generalized interaction 
operator must be of the form (36), and correspondingly the T matrix is given 
by (37). To have a bound-state corresponding to the deuteron binding energy 
Ed = 2.225MeV, in the case of the np scattering, the following condition must 
be satisfied: 

- 62 + + blMi-Ed) = 0. (79) 

Taking into account (79) and fitting the NN phase shifts in the range 0-350 
MeV, we have obtained the parameters of the interaction operator given 
by (36). The results for the S wave phase shifts are shown in Figs. (1-3). 
Values for the corresponding constants are given in Table 1. It can be seen 
from Figs. (1-3) that our model yields nucleon-nucleon phase shifts in good 
agreement with experiment despite the fact that the interaction operator given 
by (36) is rank-one separable. The phase shifts in ^Sq channels change sign as 
well as in the case of Tabakin-type potentials. At the same time, our model 
has more freedom in fitting data than the separable-potential model. 

In contrast with other phenomenological models, the form of the 



retardation effects in describing tlie A^A^ interaction. Tlie nonlocal interaction 
operator given by (36) describes the history of the two-nucleon system in 
infinitesimal time intervals during which the quark and gluon degrees of freedom 
can manifest themselves. Being exactly solvable, our model can be used for 
investigating some effects of quark-gluon retardation on low-energy nucleon 
dynamics. As we have shown, there is the correspondence between the form of 
the generalized interaction operator and the large-momentum behavior of the 
matrix elements of the evolution operator. In the nonlocal case, these matrix 
elements as functions of momenta do not tend to zero at infinity so rapidly as 
it is required by ordinary quantum mechanics, and within the Hamiltonian 
formalism this leads to the ultraviolet divergences. Correspondingly in 
the nonlocal case the two-nucleon amplitudes < p2|r(z)|pi > also do not 
decrease when |pj|^oo as it is required by Hamiltonian dynamics. Such 
a large-momentum behavior of the two-nucleon amplitudes is a consequence 
of the fact that in the nonlocal case the operator ii^non(^25 ^i) is not zero 
and must satisfy Eq.(23). Another consequence of nonlocality in time of the 
NN interaction is that for fixed momenta pi and p2 the matrix elements 
< p2|r(z)|pi > tend to zero as z ^ — oo, while, in the local case, they tend 
to < P2|^|pi > in this hmit. To illustrate this, we present in Fig. 4 the 
off-shell behavior < p2|r(z)|pi > in the limit z — oo. Thus, nonlocality in 
time of the A^A^ interaction caused by quark-gluon retardation effects gives 
rise to an anomalous off-shell behavior of the two-nucleon amplitudes. The 
off-shell properties of the amplitudes for the ordinary interaction operator and 
the operator containing the nonlocal term are qualitatively different. This is 
true even when the two interaction operators have approximately the same 
phase shifts. Such a large variation in the off-shell behavior of the amplitudes, 
even when the interaction operators are identical on-shell, can have significant 
effects on three- and many-body results [33] . This gives reason to expect that 
the anomalous off-shell behavior of the two-nucleon amplitudes can also have 
significant effects on nucleon matter properties. From this in turn it follows 
that the quark and gluon degrees of freedom can play a significant role in 
low-energy nucleon dynamics. 



VII. Conclusion 

As we have shown, the GQD provides the extension of Hamiltonian dynamics 
which can describe the evolution of a quantum system with confined degrees 
of freedom that cannot be associated with some observables. We have shown 



deal with two infinitely separated momentum (or length) scales, one of which 
is relevant for the observable degrees of freedom, and another is relevant for 
unobservable ones. This gives rise to a peculiar dynamical situation: In the 
large-momentum limit the matrix elements of the evolution operator do not 
decrease so rapidly as it required by ordinary quantum mechanics, and this 
results in the lack of continuity of the evolution operator, and correspondingly 
in the fact that the evolution of the system is not governed by the Schrodinger 
equation. The advantage of the GQD consists in the fact it permits the 
solution of the evolution problem in this case. As we have shown, this open 
new possibilities for describing low-energy hadron dynamics. This dynamics 
can be described in terms of the hadronic degrees of freedom, while the 
generalized interaction operator accounts for the coupling between the space 
describing hadron states and the space describing the states, in which the 
quark structure of hadrons manifests itself. This operator can be extracted 
from quark-gluon dynamics, or can be constructed by some quark models. 
This has been illustrated by the example of the dynamics of nucleons with the 
internal structure described by a constituent quark model. This dynamics has 
been shown to be governed by the generalized equation of motion (9) with 
nonlocal interaction operator Hint{t2^ti) given by (75). 

We have constructed a model which is an extension of the ordinary 
separable-potential model to the case where the A^A^ interaction is nonlocal in 
time, and the form of the interaction operator takes into account retardation 
effects from quark confinement. Despite its simplicity the model yields the 
NN phase shifts in good agreement with experiment. At the same time, the 
nonlocal interaction operator constructed within the model can be used for 
describing the SR part of the NN interaction. 

The main conclusion that may be drawn from the present work is that 
quark-gluon retardation can have significant effects on low-energy hadron 
dynamics We have shown that one of the important consequences of quark- 
gluon retardation effects is the anomalous off-shell behavior of the two-nucleon 
amplitudes: The elements of the two-nucleon T matrix as functions of momenta 
have the large-momentum behavior which leads to the ultraviolet divergences 
in the Hamiltonian formalism. Since the off-shell properties of two-particle 
amplitudes are crucial for the three- and many-body results, the anomalous 
off-shell behavior of the two-nucleon amplitudes can have substantial effects 
on the dynamics of many nucleon systems, and on the properties of nuclear 
matter. Because such an anomalous behavior of the two-nucleon amplitudes 
takes place only in the case, when the interaction operator (22) contains 
the nonlocal term H^^iU.tA that account for the couoline with the soace 



describing states, in which the quark structure of hadrons manifests itself, the 
above means that one cannot ignore this short-ranged term in constructing 
effective operators of the NN interaction. This is especially important for 
solving the many-nucleon problem. 

Acknowledgments 

We would like to thank W. Greiner, I. N. Mishustin and W. Scheid 
for helpful discussions and valuable comments. R.Kh.G. would like 
to acknowledge the hospitality of Institut fiir Theoretische Physik der 
Justus-Leibig-Universitat,Giessen, where part of this work was completed. 



References 

[1] S. Pal, M. Hanauske, I. Zakout, H. Stocker, and W. Greiner, Phys. Rev. 
C 60, 015802 (1999). 

[2] A. Faessler, F. Fernandez, G. Liibeck and K. Shimizu, Phys. Lett. B 
112, 201 (1982); Nucl. Phys. A402, 555 (1983). 

[3] K. Brauer, Amand Faessler, F. Fernandez and K. Shimizu, Z. Phys. A 
320, 609 (1985). 

[4] M. Oka and K. Yazaki, Phys. Lett. B 90 (1980) 41; Nucl. Phys. A402, 
477 (1983). 

[5] D.R. Entem, A.I. Machvariani, A. Valcarce, A.J. Buchmann, Amand 
Faessler and F. Fernandez, Nucl. Phys. A602, 308 (1996). 

[6] V.I. Kukulin, V.N. Pomerantsev, and A. Faessler, Phys. Rev. C 59 3021 
(1999). 

[7] Yu.A. Simonov, Nucl. Phys. A463 231 (1984). 

[8] Fl. Stancu, S. Pepin, and Ya. Glosman, Phys. Rev. C 56, 2779 (1997). 

[9] D. Bartz and Fl. Stancu, Phys. Rev. C 59 1756 (1999). 

[10] C.-F. Qiao, H.-W. Huang, and K.-T. Chao, Phys. Rev. D 54, 2273 
(1996); Phys. Rev. D 60 094004 (1999). 

[11] R. Machleidt, K. Holinde and Ch. Elster, Phys. Rep. 149 No.l, 1 (1987). 

[12] A. D. Lahiffand I. R. Afnan, Phys. Rev. C 56, 2387 (1997). 

[13] Yu.S. Kalashnikova, I.M.Narodetsky, and V.P.Yurov, Yad. Fiz., 49, 632 
(1989). 

[14] Yu.A.Simonov, Phys. Lett. B 107, 1 (1981). 

[15] A.G. Baryshnikov, L.D.Blokhintsev, I.M.Narodetsky, and D. A. Savin, Yad. 
Fiz. 48, 1273 (1988). 

[16] A.N. Safronov, Teor. Mat. Fiz. 89, 420 (1991); Yad. Fiz., 57, 208 (1994). 

[17] Yu.A. Kuperin, K.A. Makarov, and S.P.Merkuriev, Teor. Mat. Fiz., 75, 431 
(1988); 76, 242 (1989). 



[19] M. Orlowski, Helv.Phys.Acta. 56, 1053 (1983). 

[20] B.O. Kerbikov, Yad. Fiz. 41, 725 (1985); Teor. Mat. Fiz. 65, 379 (1985). 
[21] R.Kh. Gainutdinov, J. Phys. A: Math. Gen. 32, 5657 (1999). 
[22] R.P. Feynman, Rev. Mod. Phys. 20, 367 (1948). 

[23] R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals, 
(McGraw-Hill, New York, 1965). 

[24] R. Kh. Gainutdinov and A.A.Mutygulhna, Yad. Fiz. 60, 938 (1997) 
[Physics of Atomic Nuclei, 60, 841 (1997)]. 

[25] R. Kh. Gainutdinov and A.A.Mutygullina, Yad. Fiz. 62, 2061 (1999) 
[Physics of Atomic Nuclei, 62, 1905 (1999)]. 

[26] R.F.Streater and A.S.Wightman, PCT, Spin and Statistics, And All That, 
(W.A.Benjamin, New York, 1964). 

[27] Y.Yamaguchi, Phys. Rev. 95, 1635 (1954). 

[28] F.Tabakin, Phys. Rev. 174, 1208 (1968). 

[29] S.K.Adhikari and L.Tomio, Ann. Phys. (N.Y.) 235, 103 (1994). 

[30] G.Rupp and J.A.Tjon, Phys.Rev. C 37, 1729 (1988). 

[31] S.K.Adhikari and L.Tomio, Phys.Rev. C 51, 70 (1995). 

[32] V. G. J. Stoks, R. A. M. Klomp, M. G. M. Rentmeester, and J. J. de 
Swart, Phys. Rev. G 48, 792 (1993). 

[33] R. Machleidt, F. Sammarruca, and Y. Song, Phys. Rev. G 53, 1483 
(1996). 



TABLE 1. The parameters of the interaction operator obtained 
by fitting the NN date, p = IMeV'^. 



partial wave 


a 


A 


b ■ p 


d ■ p 


Ci- p^ ^ 


62 • p^"^" 


{np) 


0.499 


200.2 


433.8 


766.2 


0.015 


7.538 • 10-4 




0.499 


406.8 


356.3 


3.651 • 10*^ 


3.086 


1.779-10-8 




0.499 


134.7 


371.7 


6.763 • 10^ 


0.421 


9.564 • 10-^ 



Figure captions 



Fig.l. Phase shifts (soUd hne) in the chanel for np scattering, 
compared to the experimental data (points) [32] . 

Fig. 2. Phase shifts (solid line) in the ^5*0 chanel for np scattering, 
compared to the experimental data (points) [32] . 

Fig. 3. Phase shifts (solid line) in the chanel for pp scattering, 
compared to the experimental data (points) [32] . 

Fig. 4. The off-shell behavior of < p|r(2;)|p > in the chanel for 
np scattering. The solid curves corresponds to the model with generalized 
interaction operator (49), compared to the model with Yamaguchi potential 
with parameters given in [30] (dashed line). 



